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Abstract 

We present a practical and accurate density functional for the exchange-correlation energy of 
electrons in two dimensions. The exchange part is based on a recent two-dimensional generalized- 
gradient approximation derived by considering the limits of small and large density gradients. 
The fully local correlation part is constructed following the Colle-Salvetti scheme and a Gaussian 
approximation for the pair density. The combination of these expressions is shown to provide 
an efficient density functional to calculate the total energies of two-dimensional electron systems 
such as semiconductor quantum dots. Excellent performance of the functional with respect to 
numerically exact reference data for quantum dots is demonstrated. 
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1. INTRODUCTION 



The practical capability of density-functional theory [l| (DFT) to capture many-particle 
properties of physically and chemically relevant systems crucially depends on the availability 
of good approximations for the exchange-correlation energy functional. Significant advances 
have been achieved by means of, e.g., local and semi- local approximations, orbital function- 
al, and hybrid functionals [2] . Most of these constructions have been carried out considering 
three spatial dimensions (3D) explicitly, which is natural in view of atoms, molecules, and 
solids. Consequently, however, the field of low- dimensional physics has been left, for the most 



s developed particularly 

#Je|. 



part, outside the coverage of DFT. For example, density functiona 
for 3D fail when applied to (quasi-) two- dimensional (2D) systems 

At present, 2D nanoscale structures have a rich variety including, e.g., quantum Hall bars, 
semiconductor quantum dots and rings, and metal-oxide-semiconductor-based field-effect 
transistors 3). It should be noted that the building block here is the quasi-2D electron 
gas, which can be treated, however, by a "pure" 2D approach, i.e., on a 2D grid with 
Coulombic electron-electron interactions. More explicitly, the degrees of freedom along the 
third dimension have been suppressed and the system is effectively 2D. Then, the influence 
of the surrounding host material is taken into account with the effective-mass approximation 
manifesting itself as an effective mass and a dielectric constant in the 2D Hamiltonian. 
Within DFT, 2D systems are usually dealt with the 2D local-density approximation 



(LDA), which is accurate in terms of total energies in many cases 10 



111 ] . Problems are 



expected to arise in highly inhomogeneous systems, or when considering the very strong 
interaction regime [12], or close to dimensional crossovers [l^. Developments in 2D beyond 
the LDA 
exchange 



rave been carried out recently, leading to accurate density functionals for both 
If]] and correlation llS ] - Good results have been obtained also from orbital 
functionals designed to follow the dimensional crossover [6|, simple orbital-free total-energy 
functionals 19], as well as from DFT specialized for strongly interacting electrons 12]. 
In this work we present a 2D density functional combining the 2D generalized-gradient 



approximation [15| (GGA) for the exchange and a local density functional for the correla- 
tion [18]. The functional is shown to be accurate not only for the exchange and correlation 
energies, respectively, but also for the total energy. Thus, it could be readily used for 
the calculation of, e.g., chemical potentials and addition energies of quantum-dot systems 



2 



- quantities directly available from Coulomb-blockade transport experiments [20|. As a 
numerically efficient method, the functional is also a promising tool for numerical studies 
on confined quantum Hall systems typically dealt with single-particle theories or with the 



Thomas-Fermi approximation 



2l|. 



2. THEORY 



We focus here on the Kohn-Sham (KS) scheme within DFT to obtain the ground-state 
total energies E tot and densities n(r) of a 2D system containing N Coulomb-interacting 
electrons. The total energy can be written, in Hartree atomic units, as a density functional 

E[n] = T s [n] + E H [n] + E xc [n] + J d 2 r v(r)n(r), (2.1) 

where T s [n] is the KS kinetic energy functional, -EhH is the Hartree energy, v(r) is the 
external scalar potential, and -Ex C [ra] is the exchange-correlation energy functional. The 
latter can be decomposed into the exchange and correlation parts as S xc [n] = E x [n] + E c [n]. 
In the following sections we consider approximations for these parts separately. 

It should be noted that for E x [n] we use an approximation which depends on the spin den- 
sity and thus also applies to spin-polarized states within the spin-DFT formalism, whereas 
for E c [n] the approximation depends only the total density, and its derivation suggests that 
is suited for spin-unpolarized systems. This restriction for the correlation calls for a spin- 
dependent extension. 



2.1. Exchange 



The following 2D-GGA is based on the idea of Becke [22|, [23j for 3D systems, which was 
extended to 2D in Ref. [15|. The spin-density (n a ) functional for the exchange energy can 
be expressed as 

E x [n a ] = -7T Y, I dV Mr) [ds h K>a (r, s), (2.2) 

<T=U 

where h x>(T (r, s) is the cylindrical average with respect to s = r' - r (with s = \r' — r|) of the 
exchange- hole (or Fermi-hole) function h XjC7 (r, r') around the reference point r. Note that 
here we use a definition of the exchange energy for which the corresponding hole function 



3 



is positively defined. The Taylor expansion of h Xj(T (r, s) with respect to s up to the second 
order gives 



h, 



n„ r 



C°(r)s 2 + ... 



where 



V 2 n„(r) 



2t „ w + ^ 

2 



+ 2 



(2.3) 
(2.4) 



25|. 



is the local curvature of the exchange hole around the given reference point r [fj], Q, |24 
Here r CT is twice the kinetic-energy density and j PjCr is the paramagnetic current density. 

In the small density-gradient limit (SGL), we take the homogeneous 2D electron gas 
(2DEG) as the reference system. When the inhomogeneity is small, we may write 



, . { [1 + a a (r)s 2 + 6 CT (r) S 4 + . . .] hff G (s), if k F , a s < 



if /cfo-s > z. 



(2.5) 



Here ^° EG is the exact exchange-hole function for the 2DEG [26], kp^ = y/^nn^ is the 
2D Fermi momentum, and where z is the first zero of J\, i.e., the Bessel function of the 
first kind appearing in hff G . Comparing Eqs. (Q and fl23|) . using the 2D Thomas-Fermi 
expression for r rT) and enforcing the exact normalization of h x>a lead to expressions for a a 
and b a (see Ref. |l5| for details), and finally the SGL can be written as 



E 



SGL 



where 



1(0)7(3) - 1(1)1(2) 
48v^F [ 1(3) 

"dyy m J!(y). 



d 2 r 



|Vn CT (r)| 2 



I(m) 



(2.6) 



(2.7) 



In the large density-gradient limit (LGL) the density gradient dominates over the other 
terms in h Xj(T , and, secondly, the exchange hole vanishes at large s following a decay function 
denoted here as F. Thus, we may write 



1 \VnJr) 



F(a a (r) s) 



(2.8) 



where F(y) = e~ y2 corresponds to a Gaussian approximation for the exchange hole. The 
parameter a a can be determined by enforcing again the exact normalization constraint for 
h x ,v The resulting expression for the LGL of the exchange energy in Eq. (12. 2 ft is 



E 



LGL 



2 3/2 



n, 



3/4 



r) 



(2.9) 



where 



G(m) = f°°dy y m e~ y2 . (2.10) 
Jo 

Finally we interpolate the results for the SGL [Eq. (12. 6p ] and LGL [Eq. (12. 9p ] leading to 
the 2D-GGA expression for the exchange energy, 



<r=U 1 



nj 2 (v) 



1 + 7 



|Vrt CT (r) 



n|(r) 



3/4 



where 



LDA 



3V? 



d 2 r n^. /2 (i 



(2.11) 



(2.12) 



is the 2D-LDA exchange energy. The parameters ( and 7 are fitted to a set of parabolic 
quantum dots yielding ( = 0.003317 and 7 = 0.008323 [15(. Interestingly, these values are 



close to the parameters found by Becke by fitting a series of noble-gas atoms 



23]. 



2.2. Correlation 



Here we review the derivation of the correlation energy functional presented in Ref. 18 ]. 



The correlation energy can be expressed as 



18 



27 



28] 



E c = /dV Jdh „ 2 , SD (r, s f^>~^> , (2.13) 

where n 2i so( r , s ) is the pair density calculated using a single Slater determinant (SD) gener- 
ated within DFT from the occupied KS orbitals [see Eqs. (12.151) and (13.11) below]. Further, 
the quantity 

(p(r, s) = [1 - $(r)(l + as)} e^ 2 ^ 2 (2.14) 

describes the correlation between electron pairs. Here s = ri — r 2 and r = (r x + r 2 )/2 
represent the relative and center-of-mass coordinates of a representative electron pair, re- 
spectively. The quantities a, (3, and $ act as correlation factors to be determined as follows. 

• First, a can be found by considering the cusp condition for a 2D singlet many-body 
wavefunction. This corresponds to a situation when two electrons are brought to the 
same point. Application of the exact result of Rajagopal et al. [29( to the model wave 
function 

*(ri<7i, r N a N ) = *sD(riO"i, r N a N ) JJ [1 - (p{r h r^)] (2.15) 

n 

- from which the correlation energy given in Eq. (I2.13P is obtained [18j - yields a — 1. 



Second, we introduce /3 as a local quantity determining the local correlation length. We 
can estimate this length by comparing the "correlation area", i.e., area integral over 
the exponential part of Eq. ( I2.14p . to the area of a circle enclosing, on the average, one 



electron, i.e., 7rr 2 , where r s (r) = l/y7rn(r). The comparison yields /3(r) = qyn(r), 
where we introduce q as a fitting parameter. 

Third, the integral of (p(r, s) [Eq. f)2.14p ] over the relative coordinate s must vanish 



fsee Refs. 



28 



30| for details). This leads to a relation $(r) = /?(r)/ \p(r) + y/n/2]. 



To further simpli 
sian approximation 



30 



he expression for the correlation energy [Eq. f)2.13p ]. we use a Gaus- 



31| for the SD pair density, 

n 2 , S D(r,s)=n 2iSD (r)e-^ 2 W. (2.16) 



Applying the exact sum rule on the pair density yields 

«SD(r) = JfZ\ j^ 3 n 2,SD(r,s) = ^^n 2i s D (r)7 2 (r) , (2.17) 

where we can use the well-known relation in the SD case: n 2t sr>( r ) — \ n sD( r )- Further, 
we can associate the SD density ti-sd with the density in the second item of the list above. 
Taken together, and performing the integration over s in Eq. (12. 13ft . leads to the correlation 
energy 

E c [n] = Jd 2 rn{v) e c (r) , (2.18) 
where e c (r) is the local correlation energy per electron having an expression 

= *sv*m mr) _ !, + mm^A 

V ; 2g 2 \2 v / 2T^ L K J J 2 + c 

+ 4£M [0(r) _ 1] + iW\ (2 . 19) 



4/3(r)(2 + c) 3 / 2 \Zl~-j-~c J l + c_ 

with c = it / [2(N — l)q 2 ]. We point out that this expression includes an ad-hoc modification 
in the first term, [$(r) — l] 2 — > [$(r) — 1]. This modification is introduced to better repro- 
duce the reference values of numerically accurate correlation energies for a set of quantum 
dots and the 2DEG Q- 

Equation (I2.18P with Eq. (I2.19P define an explicit density functional for the correlation 
energy. It is self-interaction free (in the sense that it is equal to zero for one particle systems, 
N — 1) and depends on a single fitting parameter q (see below). We point out that the 



functional is not size-consistent due to the nonlinear dependence on N. In practice, however, 
this does not affect the performance of the approximation when considering finite 2D systems 
with fixed N, e.g., semiconductor quantum dots. 

The fitting parameter q = 3.9274 is chosen to reproduce the exact correlation energy 
E c = E tot — E^ x for the singlet state of a two-electron parabolic quantum dot (2D harmonic 
oscillator) with the confining strength u> = 1. Here E to t = 3 is Taut's analytic result 32 1 
and Ef xx is the exact-exchange result. For consistency, and in order to provide a predictive 
approximation, the fitting parameter q is then kept fixed for all systems. 



3. RESULTS 

The sum of Eqs. (12. lip and (" 12 . 1 8[) is the approximation we employ for the E xc . We test 
this functional in the calculation of ground-state total energies of spin-unpolarized, closed- 
shell quantum dots. Thus %(r) = n±(r) = n(r)/2 and, as a consequence, we can restrict 
our calculations to the standard DFT scheme. The total energy is then obtained from Eq. 
(12. ip in connection with solving self-consistently the KS equation 

1. 



<Pi(r) = £i(pi(r). (3.1) 



The KS potential is given by t>Ks( r ) = vn(r) + v xc {r) + t>(r), where the Hartree potential is 
computed as 

v n (r)= /dV^l (3.2) 



|r- r, 

the exchange- correlation potential is obtained from 

5E 

«*(r) = ^r y (3.3) 

and v(r) is the given external potential. The KS orbitals <fi(r) provide the the total density 
as n(r) = YliLi \Pi( r )\ 2 {N/2 electrons for each spin-channel), and e» are the corresponding 
KS energies. 

Below, we compare the reference results available from (numerically) exact calculations 
to self-consistent DFT results obtained using the real-space octopus code 
proposed functional and with the LDA, respectively. 



33] with the 
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TABLE I: Total energies (in atomic units) for parabolic quantum dots. We compare numerically 
exact results (-E^t) *° our functional (£^ e t re ) and to the local-density approximation (El^ A ). The 
last row contains the mean percentage error A. 



N 






rphere 


^tot 


2 


1 


3* 


3.026 


3.066 


2 


1/4 


0.9324t 


0.936 


0.952 


2 


1/6 


2/3* 


0.668 


0.682 


2 


1/16 


0.3031+ 


0.300 


0.308 


6 


1/1.89 2 


7.6001* 


7.629 


7.632 


6 


1/4 


6.995+ 


7.009 


7.012 


6 


1/16 


2.528+ 


2.528 


2.534 


12 


1/1.89 2 


25.636* 


25.72 


25.67 


A 






0.42% 


1.2% 



* Analytic solution by Taut from Ref. 32. ' CI data from Ref. 34. * Diffusion QMC data from 
Ref. 35. 



3.1. Parabolic quantum dots 

First we consider a set of 2D parabolic (harmonic) quantum dots, where the external 
confining potential in Eq. (12. ip is given by v(r) = u 2 r 2 /2. Table [J shows the total energies 
for N = 2 ... 12 with various confinement strengths u. The reference data E r t f t have been 



collected from analytic results by Taut 32] , configuration-interaction calculations by Rontani 
et al. 34J, and diffusion quantum Monte Carlo (QMC) calculations by Pederiva et al. [35[ 
Overall, both LDA and our functional perform very well with respect to the reference data, 
the mean percentage errors being ~ 1.2 % and 0.42 %, respectively. However, in view of 
the excellent performance of the LDA in terms of total energies, it is remarkable that the 
present functional reduces the error further by a factor of three. Moreover, it should be 
noted that the LDA total energy has the well-known error compensation from exchange and 



correlation energies, respectively 



3. 



In the results shown in Table [U for example, the LDA 



overestimates (underestimates) the exchange (correlation) energ y by 10 ... 20 %, whereas the 
corresponding errors in our functional are significantly smaller jisl Ibsl . 
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As seen in Table HJ the total energy of the largest system with N = 12 is obtained by the 
LDA more accurately than by our functional. This raises a question whether the LDA would 
considerably outperform our functional in the important large- N limit. Therefore, we also 
tested the filling-factor v — 2 state of a N = 48 quantum dot, for which relatively accurate 
variational QMC data is available (see Fig. 8 in Ref. [36|). Remarkably, the relative errors 
of the present functional and the LDA are only ~ 0.1 % (note that in Ref. 36] a different 
DFT code was used producing a slightly larger error in the LDA). This test confirms that 
the present functional is valid also in relatively large 2D systems. This is expected in view 
of the good reproducibility of the 2DEG result imposed on both the exchange [see Eq. (12.5 
and correlation (see Fig. 1 in Ref. 



18j). 



3.2. Rectangular quantum dots 

Next we consider the total energies of rectangular quantum dots with a side-length ratio \ 
and a total area of 7r 2 (a.u.) enclosed by hard- wall boundaries. As the reference results we use 
the variational QMC data in Ref. J37], where the choice for the dot area was motivated by a 
rational-valued single-electron energy spectrum, Ey = (x «+j/x)/2, where (i, j) = 1,2,3,.... 
Moreover, the dot size corresponds to an realistic area of ~ 900 nm 2 , when using the effective- 
mass approximation for electrons in GaAs [371 ] . 

The total-energy results are shown in Table [III Again, the present functional is more 
accurate than the LDA: in this case the LDA error is reduced by a factor of 1.7. There 
is no clear tendency in the accuracies as a function of N. Precise determination of the 
^-dependence would require highly accurate numerical calculations and reference data. In 
this respect the present boundary conditions (hard walls) are problematic due to the huge 
number of grid points needed. 



4. SUMMARY 



We have considered a practical and accurate density functional for the energy of electrons 
confined in two dimensions. The exchange contribution is a semi-local generalized-gradient 
approximation, and the correlation part has been built on the Colle-Salvetti scheme with 
a Gaussian approximation for the pair density. We have verified that self-consistent appli- 



9 



TABLE II: Total energies (in atomic units) for rectangular quantum dots with a total area of 7r 2 
and side-length ratio x- We compare numerically accurate quantum Monte Carlo results (37] (-E^ot) 
to our functional (E^°) and to the local-density approximation (E^ A ). The last row contains 
the mean percentage error A. 



y 
A 


A 




phere 
- C/ tot 


pLDA 


1 


2 


3.273 


3.312 


3.357 




6 


26.97 


26.98 


27.10 




8 


46.79 


46.86 


46.99 




12 


103 34 


1 D3 1 


1 D3 2 




16 


178.50 


178.4 


178.5 


2 


2 


3.696 


3.674 


3.735 




4 


12.32 


12.36 


12.45 




6 


27.15 


27.25 


27.36 




8 


47.82 


47.69 


47.80 




12 


102.26 


102.1 


102.2 




16 


177.80 


178.0 


178.1 


3 


2 


4.375 


4.321 


4.403 




4 


12.99 


12.95 


13.08 




6 


26.69 


26.75 


26.91 




8 


46.35 


46.49 


46.67 




12 


103.46 


103.4 


103.5 




16 


177.37 


177.1 


177.3 



A 0.34% 0.57% 



cation performs very well for a variety of two-dimensional quantum-dot systems. On the 
average, the present functional was found to considerably reduce the error in total energies 
over the local-density approximation. Preliminary tests suggest also good performance for 
large electron numbers with a modest computational cost. Future developments include an 
extension of the correlation part to spin-polarized systems. 
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